Pinning quantum phase transition of photons in a hollow-core fiber 
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We show that a pinning quantum phase transition for photons could be observed in a hollow-core 
one-dimensional fiber loaded with a cold atomic gas. Utilizing the strong light confinement in the 
fiber, a range of different strongly correlated polaritonic and photonic states, corresponding to both 
strong and weak interactions can be created and probed. The key ingredient is the creation of 
a tunable effective lattice potential acting on the interacting polaritonic gas which is possible by 
slightly modulating the atomic density. We analyze the relevant phase diagram corresponding to 
the realizable Bose-Hubbard (weak) and sine-Gordon (strong) interacting regimes and conclude by 
describing the measurement process. The latter consists of mapping the stationary excitations to 
propagating light pulses whose correlations can be efficiently probed once they exit the fiber using 
available optical technologies 

PACS numbers: 42.50.-p,71.10.Pm,42.50.Gy 
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The Bose-Hubbard (BH) and sine-Gordon (sG) mod- 
els have been extremely successful in describing a range 
of quantum many body effects and especially quantum 
phase transitions (QPT) 1]. Cold atoms in optical lat- 
tices have been so far the most famous platform to im- 
plement these models, where the Mott insulator (MI) to 
Superfluid (SF) QPT for a weakly interacting gas in a 
deep lattice potential was observed 0, More recently, 
it was made possible to tune up the interactions between 
the atoms in the gas leading to the realization of the sG 
model and the Pinning QPT HQ 

Alternative platforms in the field of quantum simula- 
tions of many body effects involve ions for quantum mag- 
nets 0, and more recently photonic lattices for the un- 
derstanding of in and out of equilibrium quantum many 
body effects 0- The photon based ideas are motivated 
by significant advances in the fields of Cavity QED and 
quantum nonlinear optics [8] and have initiated a a stream 
of works in the many body properties of both closed and 
lossy cavity arrays [9( . More recently, a new direction has 
appeared in the field of SCPs where hollow-core optical 
fibers filled with cold atomic gases were considered [lOj. 
The strong light confinement and the resulting large op- 
tical nonlinearities in the single photon level predicted 
for similar systems motivated new proposals to ob- 
serve photon crystallization and photonic spin-charge 
separation [l2j|. 

We will show here for the first time that is possible 
to impose an effective lattice potential on the strongly 
interacting polaritonic gas in the fiber. This opens the 
possibilities for a large range of Hamiltonians to be sim- 
ulated with photons. As first examples will analyze the 
simulation of the sG and BH models. We will show that 
the whole phase diagram of the Mott to SF transitions for 
both models can be reproduced including a correspond- 
ing photonic "pinning transition" . We conclude with a 
discussion on the available tunability of the quantum op- 



tical parameters for the observation of the strongly cor- 
related phases. The latter is possible by releasing the 
trapped polaritons and measuring the correlation on the 
photons emitted at the other end of the fiber using avail- 
able optical technology. We note here that although our 
description focus on the hollow-core case, tapered fibers 
with cold atoms could also considered [l3j 

Model setup- The considered atomic level structure is 
shown in Fig. [T] The one dimensional cold atomic ensem- 
ble is prepared outside using standard cold atom tech- 
niques and is transferred in the fiber using techniques 
described in [T^. The atoms are initially in the ground 
state | a) and the hollow photonic crystal fiber is injected 
with a quantum coherent pulse E + from the left side 
while a pair of classical fields f2± are driving the atomic 
gas from both two sides (Fig. H{b)). The atomic A con- 
figuration consisting of states la, 6, c) comprises the typ- 
ical stationary light set-up [1, |Tl| . Setting the energy of 
atomic level a to be zero and the atomic levels b, c, and 
d to —uj q + A , —u! cc + S, and — lu cc — uj q + A p , the Hamil- 
tonian describing the four-level atoms, the quantum and 
classical fields, and their interaction is written as (in the 
interaction picture) 

H = - J n a dz{5cr cc + A a bb + A p cr dd 

+ {gba<7 b a + gdcVdc) {E + e ik " z + £Le - **« 
+<r bc (Q+e^ z + ^e- ik " z ) + h.c.}, (1) 



where cry = 



(z,t) = E ± = E ± (z,t) = 

X^afce^^-^e-^-^)*, and fi± = fl±(z,t) = 
Y^ k f e ±i(k-k c )z e -i(u k -u c )t i a .. are collective and con- 
tinuous operators describing the average of |i) (j\ over 
atoms with density n a in a small but macroscopic region 
around z, where i, j go over a, 6, c, and d. The fields 
E± couple the ground state \a) to excited state \b) with a 
strength given by g ba , and |c) to |d) with a strength g dc . 



2 




/ww\ 




(d) 



FIG. 1: (Color online) In (a) and (b) an ensemble of cold 
atoms with a 4-level structure is interacting with a pair of 
classical fields Q± which create an effective stationary Bragg 
grating. The photons carried by the input pulse E + coming 
in from the left are mapped to stationary excitations (polari- 
tons) which are trapped in the grating. The strong photon 
nonlinearity induced by the 4th level leads to the creation 
of a strongly interacting gas (grey areas in (c) and (d). By 
modulating the density of the cold atomic ensemble, an effec- 
tive lattice potential for polaritons/photons can be created- 
shown by red lines in (c) and (d). Tuning to the regime of 
weak interactions between polaritons and adjusting the lat- 
tice depth, the system undergoes a BH phase transition from 
SF (upper one) to MI (lower one) phase (c). In the opposite 
regime of strong interactions, the dynamics are described by 
the sG model where the addition of an even shallow polari- 
tonic potential a "pinning transition" for polaritons (d) could 
be observed. As soon as the desired correlated state is en- 
gineered, Q- is switched off, and the excitations propagate 
out of the fiber as correlated photons. The necessary correla- 
tions measurements to probe the phases of the system can be 
performed using standard optical technology on the photons 
exiting the fiber. 



The metastable state |c) and \b) are coupled by classi- 
cal, counter-propagating control fields fi± . The quantum 
field and classical field envelopes are slowly varying op- 
erators and k q and k c are the wavevectors corresponding 
to their central frequencies u q and ui c for E± and f2±, 
respectively [ill • 

The evolution of E± in the fiber is described by the 
Maxwell-Bloch equation as 



(d t ± vd z ) E± = -ivAu}E±+iV2Trn a (g ba Vab.± + 9dc^cd,±) 

( 2 ) 

where the slowly varying operators <j a b — Ca6,+ e + 
(Jab-e~ lkqZ and a cd = cr cd:+ e lfe « 2 + a cd _e~ ik i z are in- 
troduced, v is the light speed in the empty waveguide, 
v = ui q /kq and Aw is the difference between ui q and to c 

We assume that the atoms are initialized in the ground 
state | a) and the quantum field is a weak coherent state 
containing roughly ten photons. Following the standard 
methods for treating slow light polaritons analyzed in 
HQ, and setting the coupling constants g ba = Qdc = 9 
for simplicity, we introduce 'f-f, VP- as the forward- 



and backward-going polaritons. These are the propa- 
gating excitations and are defined as ^>± = cos 9E± — 

sin 9y / 2im a a ca with cosO — VL±/ V/f 
sin = gy/2Trn a / ys 



2ng 2 n a and 



mi 



2ng 2 n a . In the limit g^J2ixn a 3> 



il±, i.e. sin# ~ 1, the excitations are mostly in the spin- 
wave form, ^± — gy / 2im a E±/Q±, with a group veloc- 
ity given by v g ~ vQ 2 / (ng 2 n a Y Next we adiabatically 
eliminate the fast rotating term from Eq. ([2]) and set 
= ('i> + + <!'_) /2 as the stationary combinations. In 
the limit of a large optical depth [14J, the equations of 
motion for \fr read 



2m 



d 2 ^ + + 2x* t * 2 



(3) 



The effective mass is m = — Aoj/(2vv g ) — 
riD^a/(4Aou g ), the potential strength is V = Auiv g /v — 
AT iT>8v g n a I (4f2 2 ) and the interaction strength between 
polaritons is \ — A 2 Sr 1D u ff /(2A p ) where A = f2 2 /(f2 2 — 
5A /2) and E = (A p - S/2)/(A p - S). 

Adding a periodic potential- The nonlinear Schrodinger 
equation ([3]) determines the evolution of the trapped po- 
laritonic field ^ (z, t) as derived from the effective Hamil- 

tonian, H = J dz^ (-^-V 2 + Vj * + x J tfe# t * t *#. 
To add an effective polaritonic lattice, we induce a pe- 
riodic atomic density distribution by applying an exter- 
nal field such that the atoms in \a) are now given by 
n a = no+n± cos 2 (irn p hz) . We keep no 3> n\ which means 
the modulation is only a perturbation in the atomic den- 
sity and derive the new Hamiltonian which reads 



H = 



-— V 2 + V + Vi cos 2 (7m ph z) 
2m 



(4) 



where Vb = Au:v g /v — KTi-£,8v g no/ (AVi 2 ) is the cor- 
responding trapping polaritonic potential and V% = 
— KTi^8v g ni/ (AVI 2 ) the resulting imposed polaritonic 
lattice depth. We stress here the dependance of the effec- 
tive polaritonic lattice on both the slow light parameters 
(group velocity, trapping laser detuning and strength), 
and the modulated atomic density. Finally, we note that 
the atomic lattice modulation should be chosen to be 
commensurate to the number of the photons in the initial 
pulse for the pinning transition to occur [Bj]. This mean 
the modulation length will approximately fall in the mi- 
crowave regime as the numbers of trapped photons in the 
initial pulse is of the order of 10 and the fiber is a few 
cm in length. 

Reaching the relevant correlated regimes- The suc- 
cess of achieving a specific strongly correlated polari- 
tonic/photonic state will be characterized by the feasi- 
bility of tuning the Lieb-Liniger ratio of the interaction 
and kinetic energies 7, and the ratio of the depth of the 
polaritonic potential to the recoil energy Vi/Er to the 
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FIG. 2: (Color online) Plots of the Lieb-Liniger interaction 
parameter 7 as a function of the one-photon detuning A p /r 
and the Rabi frequency Q/T of the classical laser field (a) 
and the lattice depth Vi/Eb. as a function of fl/T (b). The 
parameters are taken as n a = 10 7 m _1 , n\ — 0.1n a , n p h = 
10 3 m _1 , and Fid = 0.2r, A = 5r, 5 = O.Oir with F ~ 
20MHz the atomic decay rate. 



relevant regimes [J In our system these two quantities 
read: 
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Assuming fixed atomic and photonic densities, both 
quantities can be controlled by tuning the one-photon 
detuning A p /T (shifting the fourth level \d) to/from to 
resonance) and by changing the strength of the control 
Rabi frequency fi. In Fig. [2] we plot the achievable 
regimes for 7 and Vi/Er as a function of the A p /r and 
fi/r of the classical control laser field for realistic pa- 
rameters. We assume a total atomic decay rate from the 
upper level T ~ 20MHz, approximately atomic density is 
no = 10 7 m _1 (10 5 atoms into a 1cm length fiber) and a 
photonic density of n p h = 10 3 m _1 (the input quantum 
light pulse containing roughly 10 photons). For these 
values, 7 and V\/Er can be tuned in the range to 
5 and from to 20 respectively, which allows for both 
the strong and weak interaction regimes to be realized 
with our trapped polaritonic gas. We note that the cur- 
rent state of the art in the number of atoms loaded in 
similar hollow-core fibres is smaller by roughly one or- 
der of magnitude or less. However, recent experimental 
progress in the field show that our requirements should 
be satisfied in the very near future [13]. The losses which 
will mainly occur due to spontaneous emission from the 
upper levels, can be estimated by including the corre- 
sponding terms in the Hamiltonian Eq. Q. In that 
case, the effective parameters will acquire an imaginary 
part which for the effective mass for example will read 
rn = —Auj/(2vv g ) — T 1 £,n a /(4:AoVg + 2iTvg), leading to a 

n 2 v r 

loss rate n = ■ These losses will set an upper bound 

on the timescales for the preparation of the states, and 
the probing of the established correlations. For the val- 
ues under consideration in Fig. 2 and typical slow light 



velocities v g of 100m/s[llj, these translate to lifetimes 
of hundreds of micro-seconds. The latter are within the 
reach of current optical measurement technology. 

Polaritonic /photonic pinning transition- We will now 
discuss the nature of the many body states generated 
by the addition of the effective polaritonic potential and 
show that a "pinning transition" for polaritons can be 
observed, similar to the one recently experimentally ver- 
ified for bosonic atoms in This polaritonic pinning 
transition is expected to transform continuously into the 
BH regime for a sufficiently deep effective lattices (large 
V\/Er) and small interactions 7. To analyze each rele- 
vant phase of the system, we will make use of the corre- 
sponding BH and sG models from many body physics[l|. 
We will also discuss the feasibility to access the whole 
of the relevant phase diagram for both cases, by simply 
tuning the optical parameters in our system. 

We first analyze the strong interaction regime 1 < 7 < 
5, and for a weak effective potential, Vi/Er < 3. This 
regime is clearly accessible in our photonic system as we 
show in Fig. 2, by appropriate tuning of the one photon 
detuning A p and the control laser strength fl. In this case 
the proper low-energy description of the system described 
in Eq. |lj). is given by the quantum sG model which 
reads [lj 



H 





(HVg 


f d i- 





(d z ef + (d z 0Y +ym ph cos (4K6) 



(7) 

The first two terms account for the kinetic and interac- 
tion energies of polaritons respectively and d z 9 and d z <p 
denote the fluctuations of the long-wavelength density 
and phase fields 9 and <f> [1[. The dimensionless param- 
eter K = hn p hir I ijnvg) is known to be related to 7 as 
K ~ 7i7V7 - 7 3 / 2 /(27r) for 7 < 10. 

On the other hand, if we tune the system to the weak 
interaction limit with small 7 < 1 and large V\ /E-r with 
Vx/Er 3> 1, the system is characterized in a good ap- 
proximation by the BH modelQ, 



H 



-J 



E {ft 



H.c. 



7T Yl Ui (" 4 



(8) 



where J/E R = 4 OV^r) 374 exp (-2^/E^j /yfR, 

U/E K = y27^(^/i? R ) 1/4 7- 

In Fig. 3, we show that by simply varying A p /r and 
n/T the whole phase diagram corresponding both to the 
sG and BH regimes can be accessed in our system for 
realistic values of the optical parameters. We plot the 
known phase transition lines corresponding to the sG and 
BH model occurring at Vi/Er = - 73/2/ (2~7r)-4, 

and (U/J) c = v^exp (2^/W/Et) 7/ [A-k{V x / Er) 1 / 2 ] ~ 

3.85 respectively [H-Q. In our case, these are probed by 
adjusting the detuning and the laser coupling accord- 
ingly. The pinning transition is expected to occur for 
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FIG. 3: The phase diagram for sG model (large 7 and small 
Vi/En) and BH model (small 7 and large Vi/En). Here x- 
and j/-axis give the range of tuning parameters A p and fi 
needed for two regimes. The upper region corresponds to SF 
phase, while the lower region corresponds to MI phase. The 
rest of the parameters are as in Fig. 2. 
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FIG. 4: The interaction and tunneling strength as functions 
of fi/r with A p = 50F for the weakly interacting gas in the 
BH regime. The red dot line at fi/r ~ 1.03388 corresponds 
to the Mott phase transition point (U/J) — 3.85. The rest 
of the parameters are as in Fig. 2. 

any value A p /r less than 20 when fi is increased to be 
larger than T which correspond to vanishing lattice and 
strong interaction regime 7 > 3.5. The BH Mott transi- 
tion will occur in the opposite weakly interacting regime 
and deeper lattices. In Fig. 4, for a specific value of 
the Ap/r = 50 corresponding to with 7ZI, we plot he 
interaction U and and tunneling strengths as a function 
of fi/r to further illustrate this case. We see that tran- 
sition occurs for fi/r ~ 1.03388 which corresponds to 
known critical point of {U/J) c ~ 3.85. 

We will now summarize the whole process consisting of 
preparation, evolution and detection of the states. First, 
the atoms are prepared in state \a) with a slightly mod- 
ulated atomic density and transferred in the fiber [rE] ]. 
Next, E+ is sent in from the left with a co-propagating 
control field fi+ turned on. Once E + completely enters 
the medium, fi + is slowly turned off, reversibly convert- 
ing E + into the pure atomic excitations in the usual slow 
light manner [11| and the medium is illuminated simulta- 
neously by fi+ and fi_, making E+ quasi-stationary and 
trapped. The quantities 7 and Vi/E-r. are then changed 
in time to reach the desired state. Finally one of the 



control fields fi_ is switched off, mapping the polaritons 
to photons and releasing the excitations. At this time, 
the spatial correlations are mapped into temporal cor- 
relations of outgoing photons, on which measurements 
are done using standard quantum optical techniques. By 
analyzing the resulting photonic spectra and correlation 
functions the state of system could be obtained. 

In conclusion we have shown that different strongly 
correlated states of photons could be created inside 
hollow-core fibers interacting with atomic gases with cur- 
rent or near future optical technology. The resulting 
states can be controllably tuned to reproduce sG and 
BH many body dynamics and also used to probe the cor- 
responding pinning quantum phase transition predicted 
by these models. The various correlated phases can be 
analyzed by standard optical correlations measurements 
on the light exiting the fiber. 
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